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Abstract 


High temperature polybenzimidazole membrane fuel cells are the focus of attention due to high CO tolerance and overcoming water managements. 
This paper develops a transient, one-dimensional mathematical model to predict CO tolerance, and validates it with experiments. Experimental 
results are measured at different temperatures. Fuel cell performance degradations with time are also measured under various fuel compositions. 
Transient evolutions of current density, H2 coverage, CO coverage, and ionic potential are shown during the CO poisoning process. The theoretical 
results show that hydrogen coverage decreases with time, reducing hydrogen oxidation reactions and dropping ionic potential loss. The effects of 
temperature, CO contents, and H, dilutions on fuel cell performance and the time to reach steady LG are all investigated. Predictions of fuel cell 
current density degradation also show good agreement with experimental results. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Proton exchange membrane fuel cells as alternative power 
sources are gaining considerable attention due to high energy 
density and their pollution free environment potential. Pure 
hydrogen can produce PEMFC highest performance, but is not 
a proper option due to hydrogen storage and production diffi- 
culties. PEM fuel cells working at low temperature encounter 
major problems such as low CO tolerance and water man- 
agement. Preferential oxidation and water—gas-sift reaction are 
usually adopted to bring CO content down to an infinitesimal 
level. High temperature polybenzimidazole (PBI) membrane 
fuel cells are receiving recent attention. Wainright et al. [1] 
applied an acid-doped PBI membrane for fuel cells. PBI is a 
basic polymer and exhibits high conductivity through doping 
with various acids or bases. Sulphuric acid and phosphoric acid 
are the most widely used doping matter to perform high mem- 
brane conductivity. Li et al. [2] investigates the dependence of 
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electrical conductivity on acid doping level of phosphoric acid- 
doped PBI membrane fuel cell. Acid doping level, water uptake, 
water drag coefficient and mechanical strength are all investi- 
gated. Proton conductivity of PBI and PBI composite membrane 
is influenced by temperature, acid doping level and humidity [3]. 
Li etal. [4] studies the influence of doping level and water uptake 
of PBI membrane fuel cells. Methanol permeation and proton 
conductivity of poly(N-ethylbenzimidazole) (PNEBD), poly(N- 
methylbenzimidazole) (PNMBI) and PBI are investigated by Pu 
et al. [5]. Ma et al. [6] establishes a complete set of membrane 
conductivity under various RHs, temperatures and acid doping 
levels. 

The operating temperature of acid-doped PBI membrane fuel 
cells is much higher than PFSA polymer membrane fuel cells. 
Several advantages make it superior to low temperature PEMFC. 
Firstly, operating temperature up to 200°C can tolerate up to 
3% CO in anode hydrogen fuel [7]. Dhar et al. [8,9] indicates 
that CO adsorption at platinum catalysts occurs easily at low 
temperature. Thus, PBI membrane makes it possible to feed 
reformed gas directly into the fuel cell at elevated temperature 
and simplifies the fuel cell system. Secondly, water management 
such as liquid water flooding and membrane dehydration can be 
neglected at elevated temperatures higher than the boiling point 
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Nomenclature 


contact area of Pt catalyst (cm? cm7?) 
concentration of reactant gas (mol cm~?) 
diffusion coefficient (cm? s7!) 
current density (A cm7?) 
conductivity of the Nafion phase 

keg hydrogen electro-oxidation rate 
(A cm~?) 


Se e DOC E 


constant 


keco CO electro-oxidation rate constant (A cm~?) 
n number of electrons 

P total pressure (atm) 

R universal gas constant (J mol`! K7!) 
T temperature (K) 

t time (min) 

x distance (wm) 

Superscripts 

in inlet at catalyst layer 

Subscripts 

CO carbon monoxide 


CL catalyst layer 
H2 hydrogen 


MEM membrane 

O2 oxygen 

ss steady state 

Greek 

E porosity 

0 coverage ratio on Pt catalyst site 

dh ionic potential (V) 

ds phase potential of solid phase of electrode (V) 


AC A Eu) activation energy change for hydrogen dissocia- 
tive adsorption near CO occupied sites (J mol7!) 

6(AGco) variation of free energy of CO adsorption 
between zero and full coverage (J mol In 


of liquid water. He et al. [3] indicates that relative humidity has 
less influence on PBI membrane conductivity than the Nafion 
membrane. Fuel cells working at high temperature can convert 
into a compact design with methanol reformer [10,11]. This is 
because nearly 100% conversion is made around 200°C through 
methanol steam reforming [11]. Reformed gas not only feeds 
directly into the fuel cell, but also provides a sufficient amount 
of heat for fuel cells. 

Several experiments analyze high temperature PBI mem- 
brane fuel cells, but only a few theoretical studies are available. 
Cheddie and Munroe [12] apply a one-dimensional model to 
predict PBI membrane fuel cell performance. Korsgaard et 
al. [13] use a semi empirical model to approach the exper- 
imental data of PBI membrane fuel cell polarization curves. 
However, theoretical investigations examining the fundamental 
transportation inside PBI membrane fuel cells are still lacking. 
Qualitative experiments are needed to optimize fuel cell sys- 
tem designs. This paper investigates PBI membrane fuel cell 


performance under various fuel compositions from both sim- 
ulation and experiments, including effects of temperature, CO 
content, and hydrogen concentration on cell performance and 
tss. The transportations of hydrogen coverage, CO coverage on 
platinum catalysts, and ionic potential loss across the MEA are 
shown in the simulation results. 


2. Experiments 


The experiments use a five-layered membrane electrode 
assembly with PBI membrane, Pt catalyst, and carbon papers. 
Total thickness of the MEA is 971 wm and a 45.2 cm? active area. 
The end plates were made by aluminum alloy with isolated treat- 
ment to give rigid support. The serpentine flow channels were 
manufactured on graphite plates with 1 mm x | mm channel 
sizes. 

Inlet gases are supplied via model PC-540 mass flow con- 
troller (Protec Instruments, Inc.). High power electronic load 
model 3316-04 is made by Prodigit Electronics Co., Ltd. 
The ITRI test station measures polarization curves and other 
experimental results. Fuel cells are operated without gas humid- 
ification, resulting in simpler operation processes and avoiding 
water management problems, such as flooding. Fuel cell oper- 
ating temperature is around 120-180 °C through YSC (GX-36) 
temperature controllers to maintain desired operating temper- 
ature. Temperature below water boiling temperature can cause 
phosphoric acid dissolution in liquid water and drop membrane 
conductivity. The time required for MEA activation is generally 
30-50 h at a given constant current density (0.2 A em" 71 with a 
maximum stoichiometric of À =2.5 during start-up for both the 
anode and cathode side. 

Fuel cell testing is performed under various gas mixtures, 
1-3% CO, 40-55% H2, 20% CO2, balanced with N2. Pure 
hydrogen and air are applied to anode and cathode, respectively 
until fuel cell reaches steady state. Then, a CO mixed stream 
replaces pure hydrogen. 


3. Theoretical model 


A three layer mathematical model developed in this work 
analyzes the transient CO poisoning process of high temperature 
PBI membrane fuel cells. An MEA schematic model is shown 
in Fig. 1, representing the anode catalyst layer, membrane, and 
cathode catalyst layer. Table 1 represents all the governing equa- 
tions. Gas species transport at the catalyst layers is described 
as 


N; = —Djyetp VC; (1) 


where N; represents flux of gas species including hydrogen, oxy- 
gen, and carbon monoxide, ect. the porosity at both the anode 
and the cathode catalyst layer and C; expresses concentration 
of gases. Protons are produced from hydrogen oxidation reac- 
tions at the anode catalyst layer and then transport toward the 
fuel cell cathode through the membrane. Current density at the 
membrane can be expressed as: 


i = —k, Vo (2) 
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Table 1 

Governing equations 

Variables ACL MEM CCL 

Cm lec Cry] = —V An — Ri, Cm = 0 Cm = 0 

Co, Co, = 0 Co, = 0 Zo Col —V - No, — Ro, 
Cco 2 [ecLCco] = —V-Nco — Rco Cco=0 Cco =0 

$ ket Vb — (Ru, + Rco) = 0 Lé ket Vp — Ro, = 0 


in which k, is the membrane conductivity and ¢ is the electrolytic 
phase potential [14]. Water inside the fuel cell can be treated as a 
vapor phase because PBI membrane fuel cells work higher than 
120°C. Thus, the represented hydrogen and carbon monoxide 
current densities can be reduced from [15]: 


diy : np, (ds — & — Uo) 
Se = Joke 0n, sinh ( E eoe ) (3) 
dico d ncoF (Ps — d — Uo) 
= 2akeco9 h 4 
F akeco®co sin ( ORT (4) 


where key and keco are the rate constants of hydrogen and CO 
electro-oxidation [16,17], s and Uo represent phase potential of 
solid phase of electrode and thermodynamic equilibrium poten- 
tial, respectively. The calculations of unsteady state Ou, and 6co 
can be referred to [15]. The forward rate constant of hydrogen Ee 
and CO adsorption-to-desorption rate ratios bec are expressed 
as: 


ZA Ep) Ac 
RT (: ap (o>) i 


SC 


En = kfpo - exp | 


(6) 


bro = bfcoo + exp | RT 


Membrane 


x=0 | | ~ x=L 


A KÉ duu 


Fig. 1. An MEA schematic model. 


The reaction rates of hydrogen and carbon monoxide are 
expressed as: 


diy, YH 
Ry, = £ £ 7 
H2 dz (2 - ( ) 
dico / Wo 
Rco = 8 
SR dz LI, (8) 


where y and n represent stoichiometric ratio and number of elec- 
trons, respectively. Oxygen reaction rate at the cathode catalyst 


layer is given by [15]: 
Co, (==) 
exp 
COnret RT 


(9) 


1 = , 
Ro, = ap — Oco)aio | 


in which 6co is the average value of CO coverage calculated 
from the anode. Table 2 lists the corresponding boundary con- 
ditions applied in the numerical computation. Table 3 shows 
the parameters used in this paper. The steady-state condition is 
defined as the relative error of current density reaches: 


eg (10) 


in which 7” represents current density at nth min. 


Table 2 

Boundary conditions 

Variables x=0 ACL/MEM MEM/CCL x=L 

Cn, Cm =Ci Nm =0 N/A N/A 

Co, N/A N/A No, =0 Cyst 
Cco Cco = Cl, Neo =0 N/A N/A 

d bi =0 Kn ett Ei = kn $ Kn eff $ = kn bi d =0 
Table 3 


The parameters used in the present theoretical model 


Temperature (7) 453K 

Total pressure (P) latm 

Diffusion coefficient of hydrogen 1.1028 (cm? s7!) 
in gas phase (Dy, ) 


Diffusion coefficient of oxygen in 0.1775 x (T/273.15)! 83 


gas phase (Do, ) (cm? s7!) 
Thickness of catalyst layer (ôcL) 20 (um) 
Thickness of membrane (ôMEM) 30 (um) 


Gas porosity in catalyst layer (ect) 0.4 
Tonic conductivity (ba) 0.09 (mho cm7!) 
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4. Results and discussion 


Higher temperature fuel cells are expected to have higher 
CO tolerance because faster chemical kinetics facilitate CO des- 
orption on the Pt catalyst. Further advantages such as without 
water—gas-shift reaction, preferential oxidation, and water man- 
agement can obviously simplify the fuel cell system and lower 
cost. Theoretical and experimental studies both investigated the 
transportation of phosphoric acid-doped PBI membrane fuel 
cells. Fig. 2 shows the transient evolution of hydrogen coverage 
in the anode catalyst layer with time under 55% H: containing 
1% CO. Fuel cell operates without carbon monoxide poisoning 
when time equals zero. Hydrogen currently occupies the highest 
reaction surface on platinum catalysts, gradually decreasing with 
time. Strong chemical bonding favors CO adsorption on the Pt 
catalyst. At t= 0, hydrogen coverage decreases from 0.52 to 0.37 
at steady state. Thus, cell performance decreases with decreasing 
hydrogen coverage. Fig. 3 shows opposite trends and depicts CO 
coverage across the anode catalyst layer. CO coverage is initially 
at zero and reaches 0.29 until steady state. 

Fig. 4 depicts variations of ionic potential profiles across 
MEA with time. Linear ionic potential profiles are shown since 
constant ionic conductivity is assumed and protons are not con- 
sumed across the membrane. Ionic potential shows nonlinear 
profiles in the cathode catalyst layer due to oxidant consump- 
tion by reduction reactions. Ionic potential loss decreases with 
time during the CO poisoning process, lowering cell current 
density. Fig. 5 represents the distribution of 64, under various 
hydrogen dilutions with 1% CO in the anode catalyst layer. At 
40% hydrogen, 6y, has the lowest value and obviously increases 
with hydrogen content. A larger amount of hydrogen can cover 
more catalyst reaction surfaces. Increased hydrogen content at 
higher hydrogen dilution can increase a greater amount of hydro- 
gen coverage than at lower hydrogen dilution. CO coverage 
under various hydrogen contents reveals opposite trends shown 
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Fig. 2. Transient evolution of hydrogen coverage profiles in the anode catalyst 
layer under 55% H: containing 1% CO at 0.6 V. 
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Fig. 3. Transient evolution of CO coverage profiles in anode catalyst layer under 
55% H: containing 1% CO at 0.6 V. 


in Fig. 6. At lower hydrogen content, CO has more opportunity to 
cover Pt catalysts. Hydrogen dilution by comparison, has much 
more significant effect on 64, than on 6co-inlet fuel containing 
CO can decrease cell performance with time. Fig. 7 illustrates 
time to reach steady state tss under different fuel compositions. 
CO concentration amounts range from 0.1 to 3% with 40 to 
80% hydrogen content. The figure shows that tss dramatically 
decreases with increased CO content. High CO concentrations 
greater than 1% show a much smaller drop rate of tss. Differ- 
ent hydrogen contents show similar trends. Inlet fuel containing 
more than 1% CO shows only slight difference in time to reach 
steady state. Hydrogen dilution effect at low CO concentrations 
does not show significant difference on time tss, but gradu- 
ally increases with CO amount. Increasing hydrogen amount, 
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Fig. 4. The distributions of ionic potential loss across MEA with time under 
55% H: containing 1% CO at 0.6 V. 
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Fig. 5. Hydrogen coverage profiles across the anode catalyst layer under various 
hydrogen contents with 1% CO at 0.6 V. 


especially at higher CO concentrations, can obviously increase 
the time tss. Greater amount of hydrogen facilitates hydrogen 
adsorption on the platinum catalysts. 

Fig. 8 plots CO and hydrogen concentration influence on 
cell current density. Theoretical results show that increased CO 
concentrations from 0.1 to 3% can dramatically drop cell cur- 
rent density due to greater CO adsorption on Pt catalysts. Also, 
CO desorption does not easily occur due to strong chemical 
bonding between CO and Pt, thus accumulating on the catalysts. 
Increased hydrogen content does not cause significant fuel cell 
performance improvement with 0.1% CO. However, hydrogen 
dilution effect becomes much more significant with increasing 
CO concentration. At 3% CO, a significant increase of fuel cell 
current density is obtained with increasing hydrogen content 
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Fig. 6. CO coverage profiles across the anode catalyst layer under various 
hydrogen contents with 1% CO at 0.6 V. 
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Fig. 7. Effects of CO contents on the time to reach steady state under various 
hydrogen dilutions at 0.6 V. 


from 40 to 80%. Similar hydrogen dilution trends are also shown 
on time tss. Fig. 9 shows percentage of energy loss of fuel cell 
power density under various fuel compositions. A wide range of 
fuel compositions, including CO and hydrogen content, are con- 
sidered as important design parameters of methanol reformers. 
This figure allows for easy reformer specification and convenient 
performance prediction of designing fuel cells. 

Experimental measurements of PBI membrane fuel cell 
polarization curves at temperature 120, 140, 160 and 180°C 
are depicted in Fig. 10. Experiments were all conducted under 
atmospheric pressure. Fuel and oxidant were fed into the cell 
directly without extra humidification. The measurements were 
performed at a stoichiometry ratio of 2.0 of the anode and cath- 
ode. Fuel and oxidant were applied with pure hydrogen and 
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Fig. 8. Effects of CO contents on current density under various hydrogen con- 
centrations at 0.6 V. 
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Fig. 9. Percentage of power density loss under various fuel compositions at 
0.6 V. 


air, respectively. Highest performance is seen at temperature 
180°C and significantly decreases with temperature drop. This is 
because membrane conductivity increases at higher temperature. 
Electrode kinetics also becomes faster at higher temperature. 
Membrane dehydration is not shown at elevated temperature of 
PBI membrane fuel cell. Fig. 11 plots simulation result compar- 
isons with experimental data. Fuel cell performance degradation 
with time is shown under various fuel compositions containing 
1-3% CO and 40-55% H2. Transient evolutions of cell per- 
formance are accurately predicted from simulations and also 
measured by experiments. Current density drops with time and 
reaches steady state in a few minutes. Good agreements are 
plotted in this figure. The difference between simulation and 
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Fig. 10. Experimental measurements of PBI membrane fuel cell polarization 
curves at temperature 120, 140, 160 and 180°C at 0.6 V. 
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Fig. 11. Simulation result comparisons with experimental data at 0.6 V. 


experiments is relative large for the case with 1% CO and 40% 
H2. Following assumptions may cause this discrepancy. First, 
hydrogen and CO coverage are assumed to be linear relation 
with current density distribution. Second, flow channel and gas 
diffusion layer are neglected in the theoretical model. This may 
influences the transport of fuel and oxidant. Results also show 
the effects of hydrogen dilution and CO concentration on current 
density and time tss. 


5. Conclusions 


The present study develops a transient, one-dimensional 
mathematical model to analyze PBI membrane fuel cells. Vari- 
ous fuel compositions are considered to realize the effect of fuel 
composition such as CO concentration and hydrogen dilution 
of the fuel cell. Chemical bonding is much stronger between 
CO and the Pt catalyst than hydrogen; therefore CO easily 
accumulates on the catalyst. Thus, hydrogen coverage rapidly 
declines with time. Higher CO content and hydrogen dilution 
also cause significant increase on CO coverage. Ionic poten- 
tial loss decreases with time during the CO poisoning process. 
Simulation results show good agreement with experimental 
data. The simulation from this work can accurately predict fuel 
cell performance under various fuel compositions and realize 
transient degradations of fuel cell performance, thus providing 
sufficient information for the designing reformer and fuel cell 
system. 
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